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J. Rioux,' Guido Burkard,^ and J. E. Sipe' 

^Department of Physics and Institute for Optical Sciences, 
University of Toronto, 60 St. George Street, Toronto, Ontario, Canada M5S 1A7 
^Department of Physics, University of Konstanz, D-78457 Konstanz, Germany 
(Dated: March 2,2011) 

Coherent control of optically-injected carrier distributions in single and bilayer graphene allows the injec- 
tion of electrical currents. Using a tight-binding model and Fermi's golden rule, we derive the carrier and 
photocurrent densities achieved via interference of the quantum amplitudes for two-photon absorption at a fun- 
damental frequency, (B, and one-photon absorption at the second harmonic, 2m. Strong currents are injected 
under co-circular and linear polarizations. In contrast, opposite-circular polarization yields no net current. For 
single-layer graphene, the magnitude of the current is unaffected by the rotation of linear-polarization axes, in 
contrast with the bilayer and with conventional semiconductors. The dependence of the photocurrent on the 
linear-polarization axes is a clear and measurable signature of interlayer coupling in AB-stacked multilayer 
graphene. We also find that single and bilayer graphene exhibit a strong, distinct linear-circular dichroism in 
two-photon absorption. 

PACS numbers: 73.50.Pz, 78.67.Wj, 42.65.-k 



I. INTRODUCTION 

The successful isolation of a single graphene sheet'"^ has 
sparked an intense research area around its unusual electronic 
and optical properties. Carriers in graphene obey Dirac's 
equation, resulting in an electronic energy-momentum dis- 
persion that is linear, with intersecting electron and hole 
bands. ^"^ At optical frequencies, the absorption per layer 
through a graphene stack is quantized in an amount written 
in terms of universal constants.^ 

Bilayer graphene has also garnered significant interest due 
to its quite different but equivalently interesting electronic 
properties. The carriers in clean, unbiased bilayer graphene 
obey a massive Dirac equation; their band dispersion is gap- 
less, quadratic at low energy and linear at high energy.^'** 

Both single and bilayer graphene are characterized by car- 
rier mobilities that are extremely high.' '' Their high optical 
conductivity and high carrier mobilities mean they could see 
applications as optically-controlled transport devices. 

For a level system subjected to coherent irradiation at a fun- 
damental frequency and its harmonic, the quantum inter- 
ference of one- and ^-photon absorption pathways allows the 
coherent control (CC) of the excitation process. This quantum 
interference control technique has been widely used to study 
systems ranging from molecules to bulk and quantum well 
materials.'" For crystalline materials, where initial and final 
states are described by Bloch states, an often-studied method 
is the use of fundamental and second-harmonic frequencies. 
The two equivalent pathways consist of two-photon absorp- 
tion of the fundamental and one-photon absorption of the sec- 
ond harmonic. The cross-term of the transition amplitudes 
contributes to an asymmetrical distribution of injected carriers 
through reciprocal space, yielding a nonzero current density. 
The ^-space distribution is controlled by attributes of the two 
coherent components of the light field: their polarization and 
a relative phase parameter. ' ' In graphene, the interference ef- 
fect for linearly-polarized light has been predicted to be signif- 
icantly stronger than in conventional semiconductors.'- Pho- 




FIG. 1. Crystal structure of (a) single-layer graphene and (b) bilayer 
graphene. The basis vectors ai and 32 define the unit cell, 70 and y\ 
are the intralayer and interlayer coupling strengths, and the sublat- 
tices are denoted by A (A') and B (B'). 



tocuiTent CC has been demonstrated experimentally in multi- 
layer epitaxial graphene,'^ carbon nanotubes, and graphite.'^ 

In this paper, the tight-binding model of graphene near the 
Dirac point is used to calculate the distributions of caniers op- 
tically injected by simultaneous irradiation with light at a fre- 
quency 0) and light at its second harmonic frequency 2 a). We 
find that the coherent adjustment of phase parameters and po- 
larizations yield a photocuiTent for linearly and co-circularly 
polarized light. The results for one and two layers of graphene 
are contrasted, and we argue that interlayer coupling could be 
probed in a Bemal stack of graphene sheets. 

The paper is organized as follows. The effective Hamil- 
tonians used for the calculations are presented in Sec. 11. 
One- and two-photon absorption coefficients in single-layer 
graphene, the resulting distributions of injected carriers, and 
the generated photocurrents due to interference are presented 
in Sec. Ill; these results hold as well for the low-energy ex- 
pansion of bilayer graphene. The full treatment of the bilayer 
and the effects of interlayer coupling on two-photon absorp- 
tion and photocuiTent injection are presented in Sec. IV. We 
summarize and discuss our results in Sec. V. 
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FIG. 2. (Color online) Reciprocal space and linear energy-crystal 
momentum dispersion of graphene near K. The basis vectors bi and 
b2 form the reciprocal unit cell, enclosing one K and one K' val- 
ley. The dispersion shows the initially empty conduction band c and 
occupied valence band v touching at the K point. The excitation 
scheme employs interference between two-photon absorption at CO 
(red arrows) and one-photon absorption at 2a) (blue arrows), leading 
to generation of charge and current. 



where m = J\j2v\. This is the "massive" Dirac equation, de- 
scribing the carriers near the K point for jv/^TiA:! <C 7i. The 
energy dispersion consists of a pair of gapless conduction and 
valence bands touching at the Dirac point (k = 0), with a 
quadratic dependence on crystal momentum.'*'^ 

Near the K' point, similar Hamiltonians are obtained by let- 
ting —kx in Eqs. (1-3); for the purposes of this paper, the 
two valleys are equivalent. 

III. SINGLE-LAYER GRAPHENE 

In this section, we use Eq. (1) to calculate one- and two- 
photon absorption coefficients, and the CC of chiral carriers 
in graphene. The velocity operator v — j-^kH, when written 
in the eigenstates basis, takes the form 

where k is the unit vector parallel to the direction of k and 

= £ X fe. 



II. HAMILTONIANS 

Single-layer graphene (henceforth simply graphene) is a 
one-atom-thick layer of carbon atoms arranged in two trian- 
gular sublattices {A, B}, as shown in Fig. 1(a). In the basis of 
the sublattices, the tight-binding model is expanded near the 
K point to yield the effective Hamiltonian 



//i(K + k) — > hvpa • k = Vf 





hku- 



hk_ 




(1) 



where vp is the Fermi velocity, a are the Pauli matrices, k 
is the crystal momentum in the plane of the crystal relative 
to the K point and k± ~ kxi: iky. The Fermi velocity can be 
expressed in terms of the sublattice hopping term 70. "^'^ The 
resulting band energies are linear in crystal momentum k and 
are shown in Fig. 2. 

The structure of Bernal-stacked bilayer graphene is 
sketched in Fig. 1(b); there are four atoms per unit cell, each 
contributing a p- orbital to the n bands. The corresponding 
4x4 tight-binding Hamiltonian, written in the basis {A, B', 
A', B } and expanded near the K point, is given by 
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(2) 



This tight-binding model includes 70, the intralayer coupling, 
and 7i , the hopping term between sublattices A' and B from 
the two different layers. A low-energy expansion yields 



^2(K + k) 



2m 



k^_ 
kl 



(3) 



A. Carrier injection 



We calculate the rate of change of the caiTier density due 



to an interaction Hamiltonian Hi, 



where A is 



the vector potential of the optical field, niQ is the free-electron 
mass, c is the speed of light in vacuum and e = — |e| is the 
electron charge, by performing a perturbation calculation up 
to second order Assuming a monochromatic field of fre- 
quency CO, we obtain expressions for the rate of injection of 
carrier density due to one- and two-photon absorption pro- 
cesses (using the Gaussian system of quantities and cgs units 
throughout): 



«2 = ^f'''{(o)E"*{(o)E''*{(o)E'\(o)E''{(o) 



(5) 
(6) 



where E is the electric field and superscripts a, b, c, and d indi- 
cate Cartesian components; repeated superscripts are summed 
over. Microscopic expressions for the tensors ^1 and £,2 are 
derived in the independent-particle approximation following 
Fermi's golden rule (FGR)."'"' For a two-dimensional crys- 
tal, we have 



fk 



n[^v^(co,k) 5[co„(k) - to], (7) 



where a)cv(k) = a)c(k) — a),,(k), 7iC0m(k) are the band energies, 
and ntv^(a),k) is the ^-photon transition amplitude between 
valence band v and conduction band c at wavevector k: 



n[,V(a),k)--^v„(k).E(a)), 
ng)(fi) k) = y ^Vo„(k)-E(a))v„„,(k).E(co) 



a),„.(k) + a)„„,(k) 



(8) 
(9) 
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where v,„„(k) indicate matrix elements of the velocity oper- 
ator, and [(»„,f(k) + a),„v(k)]/2 = co,„{k) - [a),.(k) + co] is the 
usual energy denominator appearing in second-order pertur- 
bation theory. 

By the symmetry of graphene and bilayer graphene, the 
tensors and ^2 have respectively one and three nonzero 
independent components in the xy plane: ^2"^'^, ^2^'', 
and £,2'^^ — ^2'^^ however, all our model Hamiltonians are 
isotropic, reducing £,2 to two independent terms: ^^-'^'^ and the 
linear-circular dichroism 8 — t,^"-' / E,^^^^ 

The electric field E(a)) in an arbitrary beam at normal inci- 
dence can be written as E(a)) = Ea,e''^<^{ta+y^e'^f«')/^/2, 
for an appropriate choice of orthonormal vectors and in 
the xy plane, a real amplitude Ea, and real phase parameters 
(po) and 5<Pa,- The injection rates of the carrier density due to 
one- and two-photon processes are given by 

ni=^rM|£»P, (10) 
n2 = ir\0)) \E^t (1 - 5sin2(5(p„)) . (11) 

Both are insensitive to rotation of the crystal axes with respect 
to the normal, but one-photon absorption is independent of po- 
larization, while two-photon absorption depends on the phase 
difference 5<Pa, between the linearly-polarized components of 
the incident light. 

For the linear response, we find from Eq. (1) that = ^1, 
with 

|i(a)) =2(7o//i«, (12) 

where Oq is the universal optical conductivity of graphene: 

2 

Co = gsgvjgfi' with gs = 2 and gy = 2 denoting spin and valley 
degeneracy, respectively.*^''^ For the two-photon process we 

find ^2"^^ = = ^^''^ = -^f »■ = I2, with 

^2{co) = &gsgvfie''vli2hco)-'. (13) 

Thus, for chiral carriers, 5 = — 1 and it follows from Eq. (11) 
that circularly-polarized light {5(pco — provides twice as 
much two-photon absorption as linearly-polarized light. 

B. Quantum interference of fundamental 
and second harmonic components 

In the presence of a two-color optical field with frequency 
components CO and 2a), there exist two transition amplitudes 
connecting the same initial and final states: fi[J/'(2a),k) re- 
sults from light at 2(0 to first order in perturbation, and 
Q.CV (g), k) results from light at (O to second order in perturba- 
tion. The cross-term of these amplitudes yields the CC term. 
Although this has no effect on the total number of carriers 
optically injected in centrosymmetric crystals,'^ it yields an 
injection term for the current density. This term has the form 

ja ^ T^f"'(^(o)E''*{co)E"{co)E''{2co)+c.c., (14) 

where fj/(fi)) is a fourth-rank current-injection tensor." The 
symmetry of graphene or bilayer graphene yields rjf^™, rjy'^^. 



and rjj-'^ — rjj^ as independent components; an isotropic 
model has 277"'" = r/f™ - t^^^^-". 20,21 -y^g introduce a dispar- 
ity parameter d = rjj'^^/ 77"" to characterize how the current 
injection due to linearly-polarized beams depends on whether 
the polarization axes are perpendicular or parallel: 

j = rif^ia) (e*{co) [E* {co) • E(2a))] 

-t/E*(a)) X [E*(«) xE(2co)]) +C.C. (15) 

From Eq. (1) we find that the nonzero components of the 
current-injection tensor are related by rif'^^-^ — rjf^^ = rj^^^ = 
— T]^"-" = ifjj, and thus d = —I. In the independent-particle 
approximation, f)j is purely real. An FGR derivation predicts 
equal conduction- and valence-band contributions, for a total 
current injection 

i)i{co) ^ gsgye\li2nco)-\ (16) 

We describe CO and 2co beams at normal incidence by the 
choice of fields E{co) E^e't'o'ta and E(2a)) = E2ae''<^<-t2a, 

where ha/ia = {^mlim+yuiiim'i'^'''''^'')/^^ describing two 
arbitrary normal-incidence beams. We find that co-circularly 
polarized beams (5<Pa — 5<p2a = if ) yield the current injec- 
tion with the largest magnitude: 

j = 2V2r7/(co)£2£2»m, (17) 

where ih = ^2m sin(A^ =F 2 ) ± y2(a cos (A(p =F 2 ) . The phase- 
difference parameter A^ = 2(fift) — (fyim controls the direction 
of the current; Q is the angle that separates the polarization 
axes of the fundamental from those of the second harmonic: 
Xffl = X2(bCOS0 +y2ft,sin0. Opposite-circular polarizations 
{~8(pio = S(fy2m = ±f ) yield no net current injection, while 
linearly -polarized beams (Scpo) = S(p2o) = 0) yield 

j = 2f7/(co)£^£2<i,sin{A()r))n, (18) 

wheren = e2(ucos(20)+e2(usin(20) and 620, =zxt2a- Here 
the angle 9 between polarization axes controls the orienta- 
tion of the current within the graphene plane, and A^ controls 
its magnitude. For co-linearly polarized and cross-polarized 
beams, the injected current is parallel with the direction of 
620, the polarization axis of the second harmonic. Conversely 
the injected current is perpendicular to 620 when the polar- 
ization axes form an angle of = |. Within this model, the 
orientation of the crystal axes has no influence on the current 
injection at normal incidence. 

Graphene seems to be the first material that has been stud- 
ied for which any value of 9 is equally effective at injecting a 
current. This is in contrast to materials such as GaAs, where 
one finds | rjy'''^ \ ^ | rjf"^^ | and thus a configuration with per- 
pendicular polarization axes results in a significantly weaker 
current.! '•20.21 

Figure 3 shows the A;-space distribution of the carrier- 
injection rate, n(k) = |n[v'(2«,k)+nS'(«,k)|2, at ©^(k) = 
2co. Field amplitudes are chosen such that the integrated injec- 
tion rates from one- and two-photon processes are balanced: 
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FIG. 3. Distribution n(k) of the carrier injection thiroughi recip- 
rocal space under irradiation by an optical field with components 
E{a)) and E{2a)) satisfying A(p = f . (a) Opposite-circular polar- 
ization (— 5(pa = 5(p2(B — if' 9 = 0). (b) Co-circular polarization 
(CT* light, 5<P(o = 5(p2(B = ±f , 6 = 0). (c,d) Linear polarization 
{8(pa, = Sq^ico = 0) with 620) = X and eoj = xcos9 +ysin6; (c) 
6=0 and (d) = | . The distribution in (a) result in no net current; 
the asymmetric distributions (b-d) result in net electrical currents in- 
jected in the graphene plane along x (b,c) or y (d). 



ill (2w) = n2{co)- Opposite-circular polarization of the beams 
yield the nonpolar distribution in Fig. 3(a) and no net current. 
In Fig. 3(b) both components of the two-color field have the 
same circular polarization a*. The carrier distribution fol- 
lows n(k) oc 1 +sin(A^± 0^), where (j)^ = tan^^ (ky/kx), re- 
sulting in the injection of the current given by Eq. (17). The 
charge distribution and current rotate with A^: clockwise for 
(J+ and counterclockwise for o , when viewed from z > 0. 

In Fig. 3(c,d) we show the A;-space distribution of the 
carrier-injection rate for linearly-polarized light; without loss 
of generality, ejm is taken along the x axis; 620) = x and em = 
xcos + y sin (J. Taking the phase-difference parameter to be 
A(p = |, we maximize both the cross-term in the fc-dependent 
carrier density, n(k) oc |sin (j)k + ie-'^'P sin{2(j)k - 20) |2, and 
the resulting current, Eq. (18). For co-linear polarization axes 
as in Fig. 3(c), the distribution is symmetric with respect to ky, 
while asymmetric and strongly enhanced towards positive k^, 
although with a node at = 0. The excess of positive-A;;^ car- 
riers gives rise to a net electric current along the x axis. As the 
polarization axis of the CO component is rotated by the angle 6, 
carriers are redistributed towards positive ky. At = ^, n(k) 
is symmetric with respect to A:, and the net current is along the 
y axis [Fig. 3(d)]. 



IV. BILAYER GRAPHENE 

The chiral Hamiltonian for bilayer graphene, Eq. (3), re- 
sults in the same carrier and current injection as in the pre- 
vious section if we replace ^1 2^i, £,2 8h(o£,2/Yi, and 
fji — ^ 8/ja)fj//7i; the velocity operator in the eigenstates basis 



{V|,V2,C1,C2} 




K(K') 



FIG. 4. (Color online) Band dispersion of bilayer graphene and 
breakdown of the transition amplitudes for two-photon absorption. 
Bands vi and V2 are valence bands and initially filled, cj and C2 are 
initially empty conduction bands. Bands V2 and ci form a gapless 
doublet touching at the K point; C2 and vj are split-off bands shifted 
by an energy yi above and below the gapless doublet, respectively. 
All bands are quadratic near K and linear at larger k. Transition am- 
plitudes appear in four variants: i) the gapless term (GLT) between 
bands V2 and C], ii) two- and iii) three-band terms involving exactly 
one split-off band (2BT and 3BT, respectively), and iv) the split-off 
term (SOT) between bands V] and C2. The notation {...} next to a vir- 
tual state indicates that the sum in Eq. (9) is restricted to m e {...}. 
Not shown are the 2BT and 3BT between bands V] and c'l . 



takes the form of Eq. (4) with vp hk/m, and thus the sym- 
metry properties of the injection tensors are unchanged. But 
such a treatment describes the carriers only near the K point 
for |v/r/ifc| <c 7i and leaves out important remote bands in the 
two-photon transition amplitude. A more accurate model is 
given by Eq. (2), which is also valid for band energies on the 
order of 71 . This 4x4 Hamiltonian introduces two additional 
bands, one above and one below the Dirac point, shifted by 
an energy 71 . More importantly, it gives the correct linear dis- 
persion for larger values of k. The band dispersion near K is 
shown in Fig. 4. 



A. Carrier injection 

Starting from Eq. (2), in this section we repeat the previous 
injection-tensor calculations for the unbiased Bernal-stacked 
graphene bilayer The one-photon carrier injection of bilayer 
graphene is obtained from Eq. (12) by replacing Gq with the 
bilayer optical conductivity a from Abergel and Fal'ko.^^ 

We break down the two-photon carrier injection into four 
distinct contributions. The first (a) comes from absorption 
by the gapless doublet [leftmost transition in Fig. 4, denoted 
GLT]. The second and third contributions arise from injec- 
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tion involving exactly one split-off band and contain either (b) 
only two-band or three-band amplitudes and no cross-term, or 
(c) cross-terms of two- and three-band amplitudes; two- and 
three-band amplitudes are denoted 2BT and 3BT in Fig. 4. 
The fourth contribution (d) comes from absorption where ini- 
tial and final states are split-off bands [SOT in Fig. 4]. The 
nonzero tensor components are obtained from the symmetry 
of the matrix elements involved for each contribution, which 
yields: 



,f"(t0)=|2«(«) 



;3|2i(«)+l2c(«)] 0(2^7(0-71) 

+ |2./(«)0(2r2a)-27i), (19a) 

^->'>'(«) = -|2«(tO) + [|2;,(«) + 3|2c(«)] 0(2^« - 71 ) 

-|2./M0(2;ia)-27i), (19b) 



where 0(x) is the Heaviside step function. The matrix el- 
ements appearing in ^2a and have the same symme- 
try as those appearing in graphene. However, the contri- 
butions involving exactly one split-off band (b,c) break the 
graphene result of 5 = — 1 . Indeed, if one defines a partial 
linear-circular dichroism 5, = ^2f^/^2i^^ for ^'^ch contribu- 
tion / <E {a,b,c,d}, it follows that 5a = —1, 5/, = j, 5c = 3 
and 5ci = —1. The total linear-circular dichroism will depend 
on the relative strength of the contributions ^2a-ii, which we 
now address. 

(2) 

In computing fi„, (a),k) for the bilayer, a difficulty arises 
since it is possible for the energy denominator inside the sum 
in Eq. (9) to become exactly zero. Take for example the top 
valence band as initial state v and the second conduction band 
as final state c. When the intermediate state m is the first con- 
duction band, there exists a value of k such that hco,„ (k) = 5 7i ; 
at this k the intermediate state lies precisely in-between the 
initial and final states. This leads to a resonance in the calcu- 
lated response functions at hco — Yi- To avoid this resonance, 
we let (Om — > C0,„ + iF/lh in Eq. (9). The linewidth F accounts 
phenomenologically for dephasing due to actual population 
of the intermediate state. Other linewidths could be added to 
describe the effects of disorder or interactions, but for values 
< F we find that their inclusion does not significantly modify 
our results. 

With ^2{(o) given in Eq. (13), we have 



^2a {(0)^^2(0)) 
^2h{CO)=^2{CO) 



2hco{2hco + 3Yiy 
{2h(o + 'Yi){2hco + 2yif' 
2y\ {2}i(ii + yx){2h(ii-yi) 



(20a) 



{2h(S)) 



2 



(2to + 27i) 



I2C (ft) 




(27iC(J + 27i) 
{2h(x) 



ViF2 



{2h(x)-2yiy +Y^- ) 

-^^{2n(o + yx){2h(o-yi) 
{2hco) 

2hco - 271 \ 



(20b) 



2h(0 + 2Yi {2h(0-2yiy+r^ 
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FIG. 5. (Color online) Quantities describing two-photon carrier in- 
jection in bilayer graphene as a function of photon energy for an 
intermediate state linewidth F/yi = 0.35. (a) The individual con- 
tributions t,2a-d icL-, plain red; b, long-dashed blue; c, short-dashed 
green; d, dotted orange) from Eq. (20). (b) The independent nonzero 
tensor components (^j"" (plain black) and (^2"' ' (dashed red) from 
Eq. (19). (c) The linear-circular dichroism 5 = ^2''^-^ / ^2^" ■ 



l2d(ft)) -|2(ft)) 



2n(o{2n(o-2Yi) 
{2n(o-Yif 

(2^10)- 271 )Vf2 



(20d) 



We graph the individual contributions ^2a-d in Fig- 5(a) and 
the resulting and components of the two-photon 
carrier-injection tensor in Fig. 5(b). In contrast to the linear 
absorption, where the limits of the bilayer response function at 
low and high photon energies gave the graphene result times a 
factor of 2, "2 in the two-photon response function this "factor 
of 2" rule does not hold at low photon energy: for Tia) ^ 71 , we 
■ ■ (ft)) 

; ahcd 



find 



9^^ft)^2traphene)(«)/27i- Further, using 



Eq. (3) instead yields ^2"(Mayer)(«) ^ 8^^«^2(graphe„e)(«)/7i; 

the discrepancy is explained since the derivation using the 4 x 
4 Hamiltonian includes important three-band terms from the 
split-off bands in the two-photon transition amplitude. At high 
photon energy ^2fbibyer)(®) ^ 2^|fg'^fph,„,)(ft)) as expected. 

Two features are apparent in the response tensor at the 
thresholds for absorption into the split-off bands. The first 
feature is a pronounced shoulder in both ^^^^ and ^^^^ at 
2h(£> = Y\ due to the onset of absorption involving one split- 
off band. The second feature is the resonance which occurs at 
2h(i) = 271. There the contributions ^2fo and ^2c approximate 
the real and imaginary parts of a complex Lorentzian function, 
and contribute to a peak in and a change of sign in ^^^^ . 

In Fig. 5(a,b) we have chosen the particular value F = 
0.3571 for the linewidth. We note that is independent 
of F, but ^2h-d are not. However, the graph of these quan- 
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FIG. 6. Diagrams of the four contributions fjia-cl to the current injec- 
tion in bilayer graphene [Eq. (22)]. 



titles changes quantitatively but not qualitatively when a dif- 
ferent finite value is chosen for F. We plot the linear-circular 
dichroism 5 in Fig. 5(c) for different values of F. It can 
be seen that F has very little effect on 5. The graphene re- 
sult (5 = — 1) is reproduced for 2fico < 7i and 2hco > 37i. 
For mid-frequencies, 71 < 2h(0 < 271, the dichroism increases 
(decreases) sharply at the first (second) split-off band thresh- 
old; 5 changes sign and has a maximum value ^ 0.5 near 
2h(0 = 1.571. 



B. Current injection 

There are four contributions to the current-injection tensor 
in bilayer graphene, outlined in Fig. 6: The first (a) comes 
from absorption by the gapless doublet. The second and third 
contributions arise from injection involving exactly one split- 
off band, with the second-order amplitude containing either 
(b) only two-band terms or (c) only three-band terms. Tran- 
sitions involving only the split-off bands make up the fourth 
contribution (d). Each individual process involves matrix ele- 
ments of varying symmetry, and they contribute differently to 
T]™ and 77f-": 

77««(«) = ifiia{(o) + [3if)ib{(o) + inic{(o)] @{2h(o - 7i) 

+ifij4co)@{2hco-2yi), (21a) 

77f "(o)) = -ifiiaico) + [ifiihi(o) + 3/fj/,(a))] &{2hco - 71 ) 
-ifij,i{(o)&{2h(o-2Yi). (21b) 




4 6 2 
Photon energy (2h(o/Yi) 

FIG. 7. (Color online) The current-injection tensor rjj in bilayer 
graphene for an intermediate state linewidth F/yi = 0.25. (a) The 
four contributions fjia-d from Eq. (22): Plain red, long-dashed blue, 
short-dashed green, and dotted orange curves correspond to compo- 
nent a, b, c, and d, respectively, (b) The tensor components r]^-*-''- 
(plain black) and J]/ '"*'^ (dashed red). 



The four contributions have dissimilar values of the 
parallel-perpendicular polarization disparity parameter: = 
— l,dt, = ^,dc = 3, and d^ = — I. Their magnitudes are given 



by 



T?/fl(«) = i7/(«) 



fjibico) = fjiico) 



2hco 



{2hco + Yi] 



{2h(0 + 3Yi 



2y\ {2h(0 + Yi){2h(0-Yi] 



{2hco) 



(2hco) 



(22a) 
(22b) 



Yi{2h(0 + Yi){2h(0-Yi) 
(2/10)) 

/ 1 2hco - 271 



y2h(0 + 2Yi {2hco~2Yi)+T^ ^ 

f\id{(0) = f\i{co) {2h(0-2Yif 
{2hco-Yi) 



(22c) 



(2to-27i)VF2 j ' 



(22d) 



with 77/(0) given in Eq. (16). The parameters f\ia-d are 
plotted in Fig. 7(a). The two independent components 
of r\i in the isotropic model, rj/"'^ and T]/"", are plot- 
ted in Fig. 7(b). In the high frequency limit, the current- 
injection tensor for the bilayer tends to '?/(biilyer)(®) ^ 
2^/(graphene)(®)- ^" frequency limit, we get 

'?/lbiilyer)(®) ^ 6;;t(JT7;(^'^^' )(co)/7i, in contrast to using 
the simple Hamiltonian of Eq. (3) which neglects three-band 
terms in the second-order amplitude and gave fj/(biilyer)('^) ^ 
8/ia)77«^'^fph^„^)(a))/7i (see the start of Sec. IV). 

For mid-frequencies, there is a sharp increase in r\j at the 
first split-off band edge at 2/ia) — Y\ and a sharp decrease at 
the second edge at 2/ia) = 271 . Two features are manifest as 
a consequence of these split-off band edges: (i) In the region 
7i < 2/10 < 271 the TJ/' ''^ component changes sign, (ii) For 
2/ia) > 271 the r\Y'-'^ component becomes very small. 

The main difference between the current injection in 
graphene and in the bilayer is the contribution of components 
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FIG. 8. (Color online) Results for the current injection in bilayer graphene with linearly-polarized (O and 2(0 beams. Left-hand side: The 
disparity parameter d = rj^^-' '' /rif describing the asymmetry between parallel and perpendicular polarization axes, with F/ yi = 0. 15 (dotted 
orange), 0.25 (short-dashed green), 0.35 (long-dashed blue), and 0.45 (plain red). Right-hand side: Polar plots of f{d,d) and g(Q,d), the 
angular distributions of the projections of j parallel and perpendicular to e2£o, as a function of the angle Q between the polarization vectors 
for d = —2, —1, —0.5, 0, and 0.5. The shaded circles represent unit amplitude and dashed lines represent negative projections. The graphene 
prediction, = — 1, is highlighted. 



fjih and fjjc, each with a vastly different value of the dispar- 
ity parameter: df, = ^ and li^ = 3 while in graphene d ^ —1. 
In Fig. 8 we plot the frequency dependence of d in bilayer 
graphene. The spectrum shows a constant —1 value from 
zero frequency until a sharp increase at the first split-off band 
edge at 2hco = y\, d rises with photon energy and eventually 
switches sign. At the second split-off band edge at 2/ja) = 
d reverses sign abruptly; for the range 2j\ < 2hco < 3/1 it 
takes on large negative values as rj^-^ remains finite but rjf™ 
becomes small [cf. Fig. 7(b)]. The value of d tends to —1 at 
higher photon energy. 

We now consider current injection in bilayer graphene un- 
der irradiation by co and 2(0 beams at normal incidence. 
Choosing the electric fields E{(o) and E{2(o) as in the pre- 
vious section, the current injection is given for co-circular po- 
larization of the beams {8(pa) = 8(p2a) = ±|) by 



j = (1 - d)V21m [77,™'] ElE2com 



(23) 



and for opposite-circular polarization {—5(pco = 8(p2(o = if ) 
by j = 0. In Eq. (23), the disparity parameter d only affects the 
magnitude of the current. In contrast, for linearly-polarized (O 
and 2(0 beams {8(pco = 8(p2co = 0) forming an angle 9 be- 
tween their polarization axes, different values of d lead to in- 
jected currents with different magnitudes but also with vastly 
dissimilar angular dependences: 

j = 2lm[rir-"]ElE2cosm{A(p) 

X fie,d)t2o,+8{e,d)ti^'\, (24) 

where f{e,d) = cos^ + d sin^ 6 and g{e,d) = 
2 (1 — c/) sin2(J. Thus, the current component that is 
parallel to 620) has a nonseparable dependence on 6 and d, 
whereas the perpendicular component always follows sin20. 



Polar plots of the functions f{9,d) and g{9,d) are shown on 
the right-hand side of Fig. 8 for = —2, —1, —0.5, 0, and 
0.5. Our result for graphene (d — ~l) yields a clover-shaped 
angular distribution: the cos 20 dependence in Eq. (18). For 
more (or less) negative values of d, the lobes around 6 = f 
and ^ become more (or less) important. At c/ = 0, these 
two lobes vanish; any current injected with perpendicular (O 
and 2(0 polarization axes is completely perpendicular to 620. 
For d > there are no nodes in the angular distribution. For 
> 1 the current parallel to 620) is stronger for perpendicular 
polarization axes compared to parallel polarization axes. By 
scanning the photon energy in the range Yi < 2/ifi) < 3Yi, the 
disparity parameter d and thus the angular dependence of the 
current injection in bilayer graphene vary significantly, in 
contrast with the current injection in single-layer graphene. 
We note in particular that the sharp changes in the value 
of d near 2h(0 w 71 and 2/ia) ~ 271 should be perceived 
experimentally by rapid transitions in the angular dependence 
of the cuiTents as the photon energy is scanned. 



V. SUMMARY AND DISCUSSION 

We have calculated the response tensors for one- and two- 
photon carrier injection and two-color current injection in 
graphene and bilayer graphene. We find a strong, frequency- 
independent linear-circular dichroism 5 — —l in the two- 
photon response of graphene; for the bilayer, 5 also equals 
— 1 when 2h(0 < 71 or 2h(0 > 371, and changes sign when 
7i < 2h(0 < 371. Using the optical CC technique, in-plane 
cuiTents are generated for co-circularly polarized and linearly- 
polarized beams. Such currents could be detected in ex- 
periments making use of contacts, or by detecting the emit- 
ted THz from the accelerated charges as was demonstrated 
by Sun et al. in epitaxially-grown multilayer graphene.'^ In 
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the bilayer, the dependence on the angle 9 between linearly- 
polarized light components at (O and 2(0 is strongly sensitive 
to the photon energy for Yi < ^fi(0 < 37i. This angular de- 
pendence is in sharp contrast to the prediction for a single 
graphene layer, and could be mapped out experimentally as 
a signature for interlayer coupling in epitaxially-grown mul- 
tilayer graphene samples, which are essentially thought of as 
uncoupled graphene layers. We have assumed that the Fermi 
energy is at the Dirac point, but for a nonzero Fermi energy 
Ep, our predictions hold for hco > \Ef\. When the Fermi en- 
ergy varies across an inhomogeneous sample, Ef should be 
taken as the largest local Fermi energy. 

Our description of the bilayer considers only the strongest 
of the interlayer coupling parameters, 71, and excludes the 
next- to-nearest coupling parameters 73 and 74. These terms 
break the isotropy of the model and introduce trigonal 
warping.^ However, we feel that excluding them in a first cal- 
culation is justified. The inclusion of 73 hardly changes the 
conductivity spectrum.^'* The same can be said of including 
the next-to-nearest neighbor hoping term in the Hamiltonian 
for graphene. 

A natural extension of the current model is to consider AB- 
stacked multilayer graphene samples. Koshino and Ando have 
shown that for n layers with n even, the Hamiltonian can be 
decoupled into n/2 bilayers (with n odd, (n — l)/2 bilayers 
and one decoupled single layer). Each bilayer pair has a re- 
duced coupling strength Am 71 where A,,, = 2 cos K,„ and K,„ is 
a wavevector in the stacking direction. The response of the 
multilayer system is the sum of the bilayer systems responses 
with detuned coupling strengths. In the limit of a high number 
of layers, the wavevector K becomes a continuous variable and 
there is a continuous spectrum of resonances, smearing out the 
response tensor and eventually modeling the response of bulk 
graphite. 

In the present paper, spin and valley degrees of freedom 
have contributed only degeneracy factors of 2 to the injection 
tensors. But spin- or valley-polarized currents could also be 
photogenerated if we first lifted the spin or valley degeneracy, 
for example in a sample subjected to a magnetic or pseudo- 
magnetic field, provided that the Fermi energy is chosen ap- 
propriately. Strain-induced pseudomagnetic fields have been 
suggested to lift the degeneracy of the K and K' points by 
lOOmeV.27-8 

We conclude by pointing out the differences between co- 
herent current control in graphene and in conventional semi- 



conductors such as GaAs. In two-color CC experiments in 
gapped semiconductors, the fields are typically chosen so that 
the semiconducting bandgap Eg lies between hco < Eg < 2h(0. 
Thus, one-photon absorption at the fundamental frequency is 
energetically forbidden. However, since the band dispersions 
of graphene and bilayer graphene are gapless, for a clean, un- 
biased sample there is a nonzero joint density of states down to 
zero frequency and one-photon absorption is always present. 
This raises an issue for two-color CC experiments where the 
usual best practice is to have balanced absorption between the 
first-order process at 2co and the second-order process at co. 
To achieve this, the fundamental beam is given most of the 
power. If this beam is absorbed in the linear regime, it can po- 
tentially flood the sample with carriers that are not taking part 
in the quantum interference. We note however that this has 
not led to difficulties in observing the coherent current control 
in multilayer epitaxial graphene.'^ 

The opening of a sufficiently large gap in the band dis- 
persion of graphene would completely eliminate one-photon 
absorption at (O. The weak spin-orbit coupling offers only 
a very small gap that has been calculated to be on the |xeV 
scale for graphene^'^'^" and bilayer graphene,^' while gaps in- 
duced by the substrate '^'■'^'■'^^ or by confinement^''*'-'^ (in cer- 
tain nanoribbon geometries, similar to the way carbon nan- 
otubes can acquire a gap) are typically tens of meV. Most in- 
terestingly, in bilayer graphene gap opening can also occur 
due to z-axis asymmetry between the two layers, which can 
be field-induced.^^ Field-induced gaps are tunable up to hun- 
dreds of meV. 

The problem of linear absorption of the fundamental could 
also be circumvented by a nonzero Fermi energy Ep, taking 
advantage of Pauli blocking to prevent one-photon absorp- 
tion at CO. With the relation h(o/2 < \Ef\ < hco, one effec- 
tively has the same condition as typical gapped semiconduc- 
tors. This offers even more tunability since the Fermi energy 
can be gate-controlled, and could lead to novel electro-optical 
devices making use of coherent current control. 
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